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We study even parity metric and matter perturbations of all angular modes in self-similar Vaidya 
space-time. We focus on the case where the background contains a naked singularity. Initial con- 
ditions are imposed describing a finite perturbation emerging from the portion of flat space-time 
preceding the matter-filled region of space-time. The most general perturbation satisfying the ini- 
tial conditions is allowed impinge upon the Cauchy horizon (CH), whereat the perturbation remains 
finite: there is no "blue-sheet" instability. However when the perturbation evolves through the CH 
and onto the second future similarity horizon of the naked singularity, divergence necessarily occurs: 
this surface is found to be unstable. The analysis is based on the study of individual modes following 
a Mellin transform of the perturbation. We present an argument that the full perturbation remains 
finite after resummation of the (possibly infinite number of) modes. 

PACS numbers: 04.20.Dw, 04.25.Nx 



1. INTRODUCTION AND SUMMARY. 



The Cosmic Censorship Hypothesis (CCH) was proposed by Penrose and states that naked singularities will not 
evolve from generic initial data (see Wald for detailed discussion). The motivation is obvious: information leaving 
the singularity is impossible to predict and all physical laws in the region affected by the naked singularity will 
break down. While the CCH seems intuitively reasonable, there are some notable counter-examplea, for instance: 
the Reissner-Nordstrom and Kerr timelike singularities, certain classes of self-similar perfect fluid and dust Q 
solutions, and the self-similar scalar field In general, self-similar spherically symmetric (SSSS) space-times are a 
rich source of counter-examples to the CCH. 

These naked singularities occur in space-times with very high degrees of symmetry, and thus perhaps the naked 
singularities are due to these unphysical restraints. We seek to move away from this symmetry by perturbing the 
space-time and examining if the singularity remains naked; that is we look for stability. 

In the Reissner-Nordstrom space-time, Penrose (5||| noticed that the Cauchy horizon (CH), that surface marking the 
boundary between regions where the singularity is and is not visible to observers, would be a surface of infinite blue 
shift. Building on this, Chandrasekhar and Hartle found that metric perturbations originating from the exterior 
have a flux that observers crossing the CH will measure as infinite. Thus the CH has a "blue sheet" instability, which 
prevents observers from crossing the CH and observing the singularity. This result has been firmly established in 
further work (see and @). This suggests the Cauchy horizon is perhaps unstable in other possible counter-examples 
to Cosmic Censorship. We will examine this mechanism of instability in SSSS space-times. 

In a previous paper , the authors have shown that all self-similar spherically symmetric space-times which admit 
a Cauchy horizon, and thus a naked singularity, are stable when we inject a scalar field; that is, the flux of the scalar 
field as measured by a timelike observer crossing the Cauchy horizon remains finite. This paper extends on that work 
by considering metric and matter perturbations of SSSS space-times. In related work, Miyamoto and Harada [Tol | 
showed that the CH in the general class of space-times mentioned above is unstable at the semi-classical level. 

For the perturbations we will use the formalism of Gerlach and Sengupta . They describe gauge invariant metric 
and matter perturbations of spherically symmetric space-times by a 2+2 split and multipole decornposition. This 
formalism has been used for a variety of backgrounds, for example Schwarzschild [T5|, timelike dust 13] and perfect 
fluid 14] . As we must specify the matter content we have begun by considering the self-similar spherically symmetric 
null dust (i.e. self-similar Vaidya) space-time. This space-time is of interest because it provides a mathematically 
simple (self-similar and spherically symmetric) example of an exact solution of Einstein's equation (that is, the metric 
can be written explicitly in terms of elementary functions) with an energy-momentum tensor obeying the dominant 
energy condition that admits a naked singularity for an easily specified set of parameters. Waugh and Lake 
showed the stability of the Cauchy horizon in Vaidya space-time at the level of the eikonal approximation. 

In the following section we define the background space-time and give the conditions for a naked singularity. We 
also introduce the Mellin transform and point out how the partial differential perturbation equations can be expressed 
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as parametrised ordinary differential equations: the perturbation is then given by the inverse MeUin transform (re- 
summation) of the solutions of the ODEs. The analysis of Sections 3-5 deals with these ODEs and their solutions 
(the 'modes'). In section 3, we describe the gauge- invariant formalism of Gerlach and Sengupta. While metric pertur- 
bations are primarily useful in modelling gravitational radiation which manifests in quadrupolar moments and above 
(see, for example, JQj), for completeness we will consider all multipole moments of the perturbation. 

In section 4 we begin specifying our initial data. For completeness we must take into account perturbations arising 
from the region of flat space preceding the matter-filled Vaidya space-time. We will denote the boundary between 
these two regions Af, and define Af a,s x = where a; < is flat space-time and a; > is Vaidya space-time; here x 
is a naturally occurring similarity variable. Our initial/boundary conditions are: i) finiteness of the perturbations on 
the axis r = in the Minkowski region, ii) finiteness on A/" as x t 0, iii) flniteness on A/" as a; | 0, and iv) continuous 
matching across x = 0. 

When we solve the perturbation equations for multipole mode number I > 2, we find a two-parameter family of 
solutions in the portion of Minkowski space-time. To satisfy initial condition (i), we must set one of these parameters 
to zero. This specifies the type of solutions which will solve initial condition (ii). The second initial condition restricts 
the acceptable range of the mode number s that arises from the Mellin transform of the perturbation equations; that 
is we can only consider s > 2^ s ^ I + m where m e N. In section 5 we solve the perturbation equations in the Vaidya 
region, and we find a four parameter family of solutions. However, to solve initial condition (iii), we must further 
restrict the range of s to s > 2, s ^ Z. To satisfy initial condition (iv), one of the parameters is fixed and thus we 
have a three parameter family of solutions. 

This class of perturbations satisfying the initial conditions is then allowed evolve onto the Cauchy horizon, and 
there it is found that the flux of the perturbation is finite. Thus the Cauchy horizon is stable under metric and 
matter perturbations (/ > 2), and does not display the blue-sheet instability shown by Reissner- Nordstrom space- 
time. Interestingly, when the perturbations are allowed evolve through the Cauchy horizon and onto the second future 
similarity horizon (see below), here the flux of the perturbations will diverge, meaning that this horizon is unstable 
under metric and matter perturbations {I > 2). 

For the dipole mode Z = 1, we reproduce the familiar result that the perturbation in the Minkowski region is pure 
gauge, that is we can always flnd a coordinate system in which the / = 1 perturbation vanishes. In the Vaidya region 
there is an Z = 1 perturbation which cannot be gauged away, and this evolves through the Vaidya space-time without 
divergence. 

The I = mode is spherically symmetric, and we can use uniqueness results: in the Minkowski space-time the I — 
perturbation generates a Schwarzschild space-time, which has to have its mass set to zero to satisfy initial condition 
1. In the Vaidya region, after a spherically symmetric perturbation we again have Vaidya space-time, merely with a 
new mass function and null coordinate. 

In Section 6, we consider the problem of resummation: does flniteness of the individual modes for an allowed 
range of the parameter s imply flniteness of the full perturbation given by the inverse Mellin transform over s of the 
modes? We present what we consider a plausibility argument that the answer to this question is 'yes': a complete 
rigorous proof is beyond the scope of the present paper. We conclude with a brief discussion in Section 7. We use the 
conventions of Wald jlj] and set c = G = 1. 



The matter model we consider is the in-falling null dust or Vaidya model. Its metric is found from the Schwarzschild 
metric by replacing the constant mass term m with a function of the advanced (ingoing) Bondi coordinate v: 



where dil? is the line element of the unit 2-sphere. By setting m linear, that is m = Xv, the space-time becomes 
self-similar; that is there is a homothetic Killing vector fleld 



2. 



SELF-SIMILAR VAIDYA SPACE-TIME. 




(2.1) 
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such that C^g^n = 2^^^. Thus we can introduce the similarity coordinate x 
self-similar Vaidya metric and matter tensors become 



u/r, and in (a;,r) coordinates the 



gfj^tydx^dx" = r'^ {-I + Xx) dx'^ + 2r (l - X + Xx'^) dxdr + x (2 

X Xx Xx^ 

tfivdx^dx" = pli_,£udx^dx'' = —dx^ + - — dxdr + - — :^dr'^ , 




x + Xx^)dr^ +r^dn^, 
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r = 0,v > 




FIG. 1; Conformal diagram for Vaidya space-time admitting a globally naked singularity. The region preceding A/" is a portion 
of Minkowski space-time, and the region to the future of A/" is filled with null dust. There are three similarity horizons at which 
the similarity coordinate x is null: a; = denoted Af, x = Xc shown dashed, and x = Xe shown as a double line. We identify 
X = Xc &s the Cauchy horizon, and will call x = x^ the second future similarity horizon (SFSH). The apparent horizon given 
by 2: = 1/ A is shown as a bold curve. 

where p = m(w)/87rr^ — X/Sirr^ is the energy density, = —d^j^v is the ingoing null direction, and Minkowski 
space-time is recovered in the limit A — *■ 0. 

The matter field is switched on at w = x = (which we will call the "threshold" and is denoted N in Figure 1), 
the region u < being flat. When the matter collapses to the origin of coordinates a singularity forms. We can 
describe the global structure of the space-time by analysing the causal nature of the similarity coordinate x. Null 
homothetic lines are zeroes of grr given in H2.3|l . that is x{\x'^ — x -I- 2), thus a; = is the null homothetic line pointing 
radially inward to the singularity at the scaling origin (its past null cone. A/"). If there are other positive real zeroes 
of \x^ — x-\-2, then these represent null homothetic lines emanating to future null infinity which meet the singularity 
in the past. The lowest of these zeroes is thus the first null geodesic to leave the singularity and escape to future null 
infinity: the Cauchy horizon. It is given by 

x^ ^{l-VT^)=x,, (2.5) 

and exists for < A < 1/8. Subsequent zeroes are additional "similarity horizons": they mark the transition of x 
from timelike to spacelike or vice versa. For the self-similar Vaidya space-time, there is one more similarity horizon, 



x^ — (l + Vl-SA) = Xe. (2.6) 
2A 

For < A < 1/8 these similarity horizons are distinct and the singularity is globally naked. In the notation of Carr 
and Gundlach [3| the causal structure is rTfSfTpSs, and is shown in Figure 1. For A = 1/8 these horizons coincide 
and the singularity is instantaneously (marginally) naked; we will not consider this case in this paper. For A > 1/8 a 
black hole forms, see Figure 2. 

The second similarity horizon is, in the purely self-similar Vaidya case, the last null geodesic to leave the singularity 
and escape to future null infinity, and thus can be called the event horizon. However, to have an asymptotically flat 
model, we can match across v = v+ > with the exterior Schwarzschild space-time; in this would not be 

the event horizon. Thus we will call x = x^ the second future similarity horizon (SFSH). 

Our use of the coordinates {x, r) is motivated by the fact that they are global, important surfaces are easily defined 
in terms of the x-coordinate, and the entire space-time (Minkowski and Vaidya) is covered by r e [0, oo). Then using 
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r = 0,?; > 



X — — oo 



t; < 



r = 




FIG. 2: Conformal diagram for Vaidya space-time with censored singularity, corresponding to the case A > 1/8. A spacehke 
singularity forms at r = for i; > 0. 

the Mellin-transform the globally defined partial differential perturbation equations reduce to ordinary differential 
equations; this occurs precisely because our space-time is self-similar. The Mellin transform is defined by 



with s € C. In coordinates (a;,r), partial derivatives with respect to r only appear in the perturbation equations in 
the form r-^ and so the Mellin transform of these equations involves only derivatives with respect to x and coefficients 
independent of r. This amounts to replacing g{x,r) with r^G{x; s), in the same way that Laplace transforms amount 
to replacing g{x,r) with e'^^G{x] s). The solution is recovered by performing an inverse Mellin transform, i.e. by 
integrating over the valid range of s. See Section 6 for a discussion of the validity of this procedure and details on the 
Mellin transform. The following analysis is carried out for what we will refer to as the 'modes', i.e. the parametrised 
(by complex s) functions G{x]s). 



Next we review the formalism of Gerlach and Sengupta jljj . We exploit the spherical symmetry of the background 
by performing a 2-1-2 split of the space-time into AP (manifold spanned by time and radial coordinates) and (unit 
2-sphere), and then decomposing the angular part of the perturbation in terms of spherical harmonics. 
First some notation: Greek indices represent coordinates on the four-dimensional space-time, capital Latin indices 
coordinates on M^, lower case Latin indices 5^, and we define covariant derivatives so: 




(2.7) 



3. 



PERTURBATION FORMALISM. 



gfll'^X = 0, gAB\C — 0: 9ab:c — 0, 



(3.1) 



and a comma defines an ordinary partial derivative. 



3. A. Angular decomposition. 



We write a non-spherical metric perturbation 



9v.v = ~gp,v + hf,i,{t,r,e,(j)), 



(3.2) 



where the over-tilde denotes background quantities; and similarly for the matter perturbation 



t 



V + AV(^,^^',0)• 



(3.3) 



From the spherical harmonics we can construct bases for vectors, 



(3.4) 
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and tensors, 

Yjab ; Zab =Y,a,b+^l{l + l)Yjab ] S(^a:b}, (3.5) 

where summation over the mode numbers I, m is impHed. Using these, we decompose the perturbation in terms of 
scalar, vector and tensor objects defined on M^, times scalar, vector and tensor bases defined on S"^. 

When we compute the linearized Einstein equations for a perturbed space-time decom pos ed in this way, we find 
they naturally decouple into two sectors, even and odd. Gundlach and Martin-Garcfa's [lj| definition is the most 
straightforward: that sector whose bases are in even powers of ej' is called even (or polar or spheroidal); that sector 
whose bases are in odd powers of ej' is called odd (or axial or toroidal) . In this paper we will consider only the even 
sector. 

We write the (even) metric perturbation as 

, _ f flAsY hAY,a \ /o f:\ 



and the (even) matter perturbation as 



l^Symm r^KY-fab + GY.a.b) 



( MabY At^y„ 



^^^^ - \ Symm r^M^Y^^b + At^Z^, ) ' ^^'^^ 

Thus the metric perturbation is defined by a symmetric two-tensor, a two-vector and two scalars, 

hAB, hA: G, K, (3.8) 

and the matter perturbation is defined by the same. 

At AS, At A, At\ At\ (3.9) 

When we write down Einstein's equations for this metric and matter tensor, we identify the left and right hand side 
coefficients of the scalar, vector and tensor bases given in H3.4|l . (|3.5|l . and these are our evolution equations for the 
perturbation. The great simplification is that these equations are in terms of two-dimensional objects Hab etc, and 
their derivatives, defined on the manifold A/^. This makes calculating the perturbation equations much easier as the 
connections used in calculating derivatives e.g. hA\B are those defined on the background. 

It is important to note that for ? = all of the basis objects given in H3.4|) . H3.5(I vanish except for Yjab, and for 
I = 1 Zab and S(^a:b) vauish. Thus some equations do not exist in these cases, and we must consider I = and / = 1 
separately from the more general I > 2 case. 



3.B. Gauge invar iance. 



Two space-times are identical if they only differ by a diffeomorphism (we take the passive view of a diffeomorphism 
as a coordinate transformation). There is a danger that if you add a "perturbation" 

g/iu = giiu + hf^^, (3.10) 

you are in fact still looking at the same space-time after undergoing a coordinate transformation, rather than after 
being perturbed in a physically meaningful way. To escape this problem we must only interest ourselves in those objects 
which do not change under an infinitesimal coordinate (gauge) transformation. These are called gauge invariants, and 
are the true measure of a physically meaningful perturbation. (More precisely these are identification gauge invariant, 
see Stewart and Walker [l^ for details). 

The vector field generating a gauge transformation has an even/odd parity decomposition as 

$,^UYdx^ +£,Y,adx\ (3.11) 
The gauge change induced on the tensor h^i, is 

h^jj ^ h^v ^^^v ^" ^^9^1^: (3-12) 
where an over-bar represents gauge transformed objects, and where 



(3.13) 
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with the covariant derivative associated with the background metric. Thus we can write down how all of the 
perturbation objects in /i^^ and At^^ transform p^ . 

Now we construct the gauge invariants, that is we take linear combinations of 1)3. 8|l . 1)3. 9|l to form objects which do 
not change under a gauge transformation. These are 

kAB = hAB ~ iPA\B + PB\a) 

k = K~2i^pA 

Tab = At^s - tAB\cP'^ - iA^Pc\B - iB^Pc\A 

Ta = MA~h''pc-r-'{i\/'^)G,A 

Ti = Ati - {p^ /r^){rH\/1),c + l)(i\/4)G 

T2 At2 - (r2i\/2)G 



(metric) (3.14a) 
> (matter) (3.14b) 



where pA — hA — ^r^G,A and — The perturbation evolution equations are then recast entirely in terms of 
the gauge invariants. We give these equations in Appendix A. 

There is an especially useful gauge choice we can make, called the Regge- Wheeler or longitudinal gauge. This consists 
of transforming to a specific coordinate system via 

hA - ^G.A ) Ydx^ + '^GY,a dx\ (3.15) 

in which Ha — G = Q ^ pa- Since we are measuring gauge invariants we are free to make this transformation, the 
benefit being that in this gauge the bare perturbations (|3.8() . H3.9|l and the gauge invariants match. 
When / = 0,1, we cannot construct a set of gauge invariant objects like those given above for the same reason 
that there are less equations in these sectors (see Appendix A): the vanishing of some or all of the bases given in 
(|3.4|) . H3.5I) . Thus for I — 1 modes we can at best construct only partially gauge invariant objects, and for Z = all 
remnants of gauge invariance are lost. 



Finally we must consider what to measure on the relevant surfaces in order to test for stability. Following 
Chandrasekhar we use the Weyl scalars to measure the flux of energy of the perturbations. For a detailed 
discussion on how the Weyl scalars relate to the Zerilh, Regge- Wheeler, Moncrief etc. scalars, see Lousto Stewart 
and Walker showed that the only Weyl scalars which are both identification gauge invariant, which is the sense 
described above, and tetrad gauge invariant (independent of the choice of null tetrad with which the Weyl scalars are 
defined), are the Petrov type-N terms. Furthermore, they are only tetrad and identification gauge invariant if the 
background is type-D or conformally flat. Since Vaidya space-time is spherically symmetric, and therefore type-D, 
this means there are two fully (tetrad and identification) gauge invariant Weyl scalars (for Z > 2), J^'o and 5"^^. 
These scalars represent pure transverse gravitational waves propagating in the radial inward (respectively outward) 
null directions of a spherically symmetric background. 
With a particular choice of null tetrad these scalars can be given as 

5-^o^^i'^PkAB(w''w''Y..ab), 5-^i^^n'^n''kAB(w*''w*''Y..ab), (3.16) 
2r^ 2r^ 

where n'^, m^ = r~^ w'^JO, </)) and m*'' are a null tetrad of the background and the * represents complex conjugation 
(see, for example, Nolan [21|). Tetrad gauge invariance of these objects must be carefully interpreted. In the type-D 
background, there is an obvious choice of null tetrad: we take £'^,n'^ to be the the principal null directions of the 
Weyl tensor, and take and its conjugate to be unit space- like vectors orthogonal to i?^, to complete the tetrad. 
Then J^'o and S'i'i are identification gauge invariant and also tetrad gauge invariant with respect to any infinitesmal 
Lorentz transformation of the tetrad, and also with respect to any finite null rotation that leaves the directions of 
£'^,h'^ fixed. However these scalars are not preserved under the finite boost-rotations 

e^' , ^ A-^n^\ fa e*"m, A > 0, G M 

under which they rescale as 

For the sake of boundary conditions however we must take this scale covariance into account, (see Beetle and Burko 
[2^ for a generalized discussion), and thus to first order our "master" function will be 

SP^i |(5^'o(5*4|^^^- (3.17) 
For / = 0, 1 the angular part of S'^o, (5^4 is zero, and thus these scalars vanish in these sectors. 
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4. MINKOWSKI REGION. 

For completeness our analysis must include the region of flat space proceeding the Vaidya region. In our global, 
self-similar coordinates (x,r), Minkowski space-time is given by 

g^^^dx^dx" = -r'^dx'^ + 2r(l - x)dxdr + x{2 - x)dr'^ + r^dVt^, t^^^dx^dx" = 0. (4.1) 

We seek to identify only that class of perturbations which are finite on the portion of the axis r = 0, u < 0, and on J\f, 
the surface marking the transition between flat space-time and Vaidya space-time, which we call the "threshold" . Af 
is the past null cone of the origin {x, r) = (0, 0). We consider the case where the perturbed space-time is still empty, 
that is the matter perturbation is zero. Thus the right hand sides of all perturbation equations is zero. 

4. A. I > 2 modes, Minkowski. 

The perturbation equations for Minkowski background are 

2v^{kjiB\c - kcA\B - kcB\A - gABkco'^ + dABkcDV^) - ^-^—^kAB 



-gAB 



(^2k,c ~- ^^f^fc) + 2{iAk,B +iBk,A +k,A\B ) = 0, 

{i^^v'' +d^^'')kcD = 0, fc,^-fc^c"' = 0, k^^^O, (4.2) 



where v"^ = r^^ jr. 

We write these equations in (a;, r) coordinates and perform a Mellin transform 12. 7|) over r. There are four unknowns, 
which, after the Mellin transform, we denote 

/ r"+MJa;) r'BAx) \ , s-i r ^ o^ 

= [ r^B.(V r^-'cl{x) j ' (^.3) 

We deflne a new variable D — B — xA, where for brevity we have dropped the subscript s. The scalars we wish to 
measure then become 

5*0 = -^r'-^D, <5*4 = ^r''-3(2yl + D), 5P-i = \6-^o6-^ . (4.4) 

We can decouple the system to get an ode in D, 

x{x - 2)D" + 2(1 + s + a; - sx)D' - {I + P + s - s^)D = 0, (4.5) 

which is hypergeometric, 

z(l - z)D"{z) + (7 - (a -f /3 -f l)z)D'{z) - af3D{z) = 0, (4.6) 

with a = l + l — s,p — —I — s,j — —1 — s, and z = x/2. There are two solutions near the axis x = — oo [23, denoted 
by the subscript A, 

^Di ^ {-zy~^-\Fi{l + I - s,3 + l;2 + 21; z~^), (4.7) 
= - H-z)^D, + {-zy+'{l + O(z-i)), (4.8) 

since a~(3 — 21 + lE'E. 

If we form the general solution 

D = di ^Di + d2 ^D2, (4.9) 
and from this solutions for A, B, C and we find to leading order near x — — oo (r — 0), 

S^OA - dy-^{l + 0{r)) + d2r-^-^{l + 0{r)). (4.10) 
Thus in order to make <^^'o,4, and hence SP^i, regular at the axis, we need to set ^2 = in the general solution for D. 
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Now we allow only the first solution for D to evolve up to the past null cone of the origin M, the threshold. When 
we do so we use the nature of the hypergeometric equation to write the acceptable solution at the regular axis as a 
linear combination of the two solutions on Af. That is to say, near x = this solution has the form 

^Di=d3^Di+d4^D2, (4.11) 

where t^i,t£'2 are two naturally arising linearly independent solutions of the hypergeometric equation near a; = 0. 
In finding these solutions, the relation between a, (3 and 7 is key so we must consider two cases: s G Z and s ^ Z. 

The more straightforward case is when s ^ Z. Then 1 — 7 ^ Z and we can take 

^Di= F{a,l3;r,z), Ti?2-2'~''i^(a-7 + l,/3-7 + l;2-7;^), (4.12) 
and l|nT|l holds with ^ 

_ r(i-7)r(a-/j + i) _ r(7)r(i-7)r(a-/3 + i) 
''^-r(i-/3)r(a-7 + i)' r(2-7)r(7-/3)r(a) ■ ^ 

Again the solutions for A, B, C and K can be recovered from these expressions for D. 

When we calculate the scalars due to jDi, near x — (and away from the singularity ar r = 0) we find 

(5*0 - ^3(1 + 0{x)) + d4x'+\l + 0{x)), d3(l + 0{x)) + dix'-^l + 0{x)). (4.14) 

These two scalars and SP-i will be finite on a; = iff 

Re{s) > 2, s^Z. (4.15) 

We can also see this with a general argument: the other singular point of the hypergeometric equation is at z = 1, or 
X = 2. This is the future null cone of the origin in Minkowski space-time. We would expect solutions to be regular 
here too, and from the transformation 

2Fi(a, 7; ^) = (1 - ^)^-"-''2i^i(7 - a, 7 - /?; 7; z) (4.16) 
we require i?e(7 — a — P) ~ Re{s — 2) > 0. 

The case s e Z is cumbersome and we will only summarize the results here. This case is again split accord- 
ing to the sign of 1 — 7 G Z, however there are no solutions for 1 — 7 < which are finite on J\f and thus we only 
present here the solutions for 1 — 7>0 (s-|-2>0). 

If 7 = 1 — m where m is a natural number, and a, /3 are different from the numbers 0, —1, —2, . . . , 1 — m, then there 
are two linearly independent solutions near a; = given by 25] 



tZ?i = zi-^F(a-7 + l,/3-7 + l;2-7;z), 

00 



where 

h* (n) = 7/1(0; + m + n) + tjj{(3 + m -|- n) — -0(1 -\- m + n) — 7/1(1 -|- n), (4.18) 

and 7/ — 7/;(°) = ^ is the Digamma function. The scalars (55*0,4 and P_i will be finite on a; = due to these solutions 
if i?e(s) > 2. 

If, however, a ov [3 is equal to one of the numbers 0, —1, —2, . . . ,1 ~ m, then the solution given above loses meaning; 
this will occur for s ~ I + m, where to is a natural number (> 1). In this case two linearly independent solutions are 

Ti?i = F(a,/?;7;z), (4.19) 

Ti?2-i^(a-7 + l,/3-7 + l;2-7;^)- (4-20) 

The scalars due to these solutions will diverge on a; = 0. Therefore we must not consider the modes s = I + m when 
summing to find the general solution. This does not challenge the generality of the result however, as all multipole 
mode numbers I (> 2) will still be counted. 

Thus we have found the largest class of perturbations that will have finite flux on the axis and on the threshold: all 
s G C such that Re{s) > 2, except for s = I + m where to is a natural number. These solutions will be matched 
across a; = into the Vaidya space-time, and then allowed evolve up to the Cauchy horizon. 
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4.B. / — 1 mode, Minkowski. 



While the resuh in this sector is well known, for the sake of completeness we briefly give the analysis in terms of 
the formalism defined above. 

In this sector not all the perturbation equations apply, and we can only define partially gauge invariant objects. 
Looking at H3.6|l . for ^ = 1 we find Yjab = —Y,a:b and thus 

h^^dx^'dx" = ... + r^{K - G)Yjabdx°-dx''. (4.21) 

If we let H = K — G, then the gauge transformation of H generated by S^^dx^ = £^AYdx^ + £Y,a dx"" is 

H - H = K - K ^ (G ^ G) = 2i/r^ - 2w^Ca- (4.22) 

Now we rename H as K , and we have effectively set G = with K now transforming as H given above. Thus 
PA = h-A, and we are left with this sensitivity to the angular part of the gauge transformation, 

kAB ^ kAB + [r\^/r'),A ] + ] |^ (4.23) 

k ~* k + 2^/r'^ + 2i\'^{^/r^),A. (4.24) 

We can, however, use this to our advantage. Following Sarbach and Tiglio [l^ . we look to transform into a coordinate 
system in which /c^'^ = 0. To do this we choose ^ such that 

[r'{^/r%A]^^ = -kA^. (4.25) 
Then we are free to make further gauge transformations, provided 

P(^A^),a]'-' = 0. (4.26) 
Thus we can reinstate the second scalar perturbation equation, fc^"^ = 0, as a gauge choice. 

Using this, we split the tensor perturbation equation into its trace and trace-free parts. Then the set of equations for 
Z = 1 is given by 

2v'^'{kAB\C - kcA\B - kcB\A) - -^kAB + 2(w^/c,_B +V^k,A +k,A\B ) - 9ABk,J^ = 0, 

2d^ij'^kcD = kj^' +2v^k,c kj^(}^ = k,A (4.27) 

where the first equation is the trace-free part of the tensor equation, and the second is its trace. 

Now we must consider what to solve the equations for. The scalars (5 5*0, 4 given before have an angular dependence 
which is zero for I — 1, and thus they vanish in this sector. The other options for true scalars to measure are fc^"^ 
and fc, however we have chosen a gauge in which the trace of kAB is zero, and thus the only scalar left to measure is 
k. We will show this scalar is pure gauge, that is there is enough residual gauge freedom to transform into a gauge in 
which fc = 0. In this gauge, all the components of kAB are also zero, and thus this perturbation sector is empty. 

For convenience we look at the perturbation equations in orthogonal coordinates, 

ds^ = -dt^ + dr^ + r^dn^, (4.28) 

and in this coordinate system we label components 

_ ( A{t,r) B{t,r) \ 

(,B(t,r) Ait,r) j' ^^■^^> 

since kAB is symmetric and trace-free. Now when we look at the perturbation equations, the trace equation gives A, 
and one of the trace-free equations gives B, in terms of k and its derivatives, 

A{t,r) = ^(2k,r+r(k,rr-k,tt)), 

B{t,r) = ^(2k,t+r\k,tt-k,rr),t)- (4.30) 

The remaining equations give 



r(nfc),^ = -4(nfc), 



(4.31) 
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which we solve as 

Dk^fit)r~^. (4.32) 

k should satisfy this equation with initial data fc(0, r) = ct(r), k{0, r) — (3{r) satisfying a, /3 e and k{t < 0, r = 0) 
finite. This implies f{t) — 0, and thus k solves the homogeneous wave equation 

□fc = 0. (4.33) 

k however is not gauge invariant. Looking at (|4.24|) . we see k ^ k — k + rj where 77 = 2£^/r'^ + 2v^r^{£^/r^),A ■ Thus 

Ok nk = Dk + nri, (4.34) 

but Dt] — (provided (|4.26|l holds), and thus k and 77 satisfy the same equation. Therefore we can choose 77 = —k, 
and thus we can always transform to a gauge in which A; = 0. In this gauge we also have kAB = 0, and thus the entire 
I — 1 perturbation is pure gauge. 



4.C. I = mode, Minkowski. 



The I — mode represents a spherically symmetric perturbation. As we are considering zero matter perturbation, 
we have the perturbed space-time is spherically symmetric vacuum, and thus Birkhoff's Theorem applies, that is the 
perturbed space-time is Schwarschild. We will recover the specifics using our perturbation formalism described above. 

For I — 0, Y,a = and thus Ka = G — Q. We cannot form gauge invariants and thus use K and Hab, which are 
fully dependent on gauge transformations £,i_idx'^ = S^Adx"^ as 

hAB hAB - {(.A\B + ^B\a), 

K K-2v^U. (4.35) 
We can transform to a gauge in which K = h"^A = by choosing (^a such that 

Ca'^ = ih^A, v^U = (4.36) 

In {t,r) coordinates, this means making a gauge transformation = and £,t.t = £,r,r — ^^a"^- Further trans- 

formations preserving K = h^j^ = must be of the form £^a = f{r)SA for any arbitrary function /(r). Thus of the 
remaining perturbation terms, A{t,r) is gauge invariant, whereas B{t,r) B{t,r) — f'{r), where 

The perturbation equations reduce to A^t = B,t = and rA,r — —A. Since _B is a function of r alone, we can choose 
/(r) to set B = 0. Thus the only perturbation term which cannot be gauged away is 

A^-. (4.38) 

r 

Renaming the constant c = 2m, and noting (1 — 2m /r) ^ « (1 -I- 2m/r) for m small (which is the case in this linear 
model), we have recovered the Schwarzschild line element. There is an intrinsic singularity on the axis r = 0, and 
thus this solution contradicts our initial data; therefore there is no Z = perturbation. 



5. VAIDYA REGION. 



In this section we describe metric and matter perturbations of the Vaidya solution given in H2.3|l . (|2.4|l . Initial 
conditions for this problem are that any perturbations coming from the region of flat space preceding the threshold 
should be finite on, and match continuously across, this surface. 
As before, we will need to split the analysis into the / = 0, ^ = 1, and / > 2 cases. 
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5. A. / > 2 modes, Vaidya. 

We consider the case where the perturbed space-time is that of a null dust, and thus has matter tensor 

t^, = (p + Sp){e^ + se^){£, + 6U), (5.1) 

where p + Sp is the energy density and + Si^ is the null vector of the perturbed space-time. As before we decompose 
these perturbation objects in terms of the spherical harmonics. Since we are measuring the gauge invariants given in 
(|3.14|) . we can examine in any gauge we choose. For convenience, we will use the Regge- Wheeler (RW) gauge given 
in H3.15|l . Thus we can write the matter gauge invariants as 

We define F as the perturbation in the null coordinate; that is we can define the null coordinate of the perturbed 
space-time a&V = v — F(a;, r)Y[B^ 0). It follows from the conservation equations that the null vector of the perturbed 
space-time is £^ -I- ^^^J, = —d^V. 

We write out the perturbation equations as given in appendix A, and as before perform a Mellin transform over r. 
There are now six unknowns, 

kAB = ( ""'^tjj^fi)^ X-^c[t) ) ' ^ = ^'~'Ks{x), F = r^G^ix), Sp = r^-Hl,{x). (5.3) 

A constraint equation defining 5p decouples and so we are left with five unknowns in the evolution system. We use 
the scalar equation ljA.ld|) to remove one of the metric variables, and so we can reduce the set of equations to a four 
dimensional first order linear system. We give this system in Appendix B. 

As before (dropping the subscript s), we define the variable D — B ~ xA, and the scalars to measure become 

Ax — 1 2 
From our system of equations we can decouple a second order ordinary differential equation for D, 

D"{x) +pix)D'{x) + q{x)D{x) = 0, (5.5) 

where 

, , s + \ 2A s — 3 - (s - 6)Aa; 

- + T3a^ + A.^-. + 2 ' (^-^^ 

(Ax - If {I + P + s^jXx - 1)) + A(x - 2 + 6A.T - Xx^) - ^(Aa; - 1)(1 + A(-2 + x{2Xx - 3))) 
" x(Aa:-l)2(A.T2-x + 2) ' ^ ^ 

The regular singular point at x = is the past null cone of the origin of coordinates, and is the surface over which 
we move from fiat space-time to Vaidya space-time (the threshold). When A < 1/8 the space-time has the structure 
given in Figure 1, and there are two surfaces which are regular singular points of the above equation: the first is the 
Cauchy horizon at x = Xc and the second is the second future similarity horizon (SFSH) x = Xe, the first and second 
zeroes of Ax^ — x -I- 2 respectively. Finally Xq = 1/A is a regular singular point describing the apparent horizon. 

We will begin the analysis by examining the threshold x = and ensuring our initial conditions arc met. As the 
analysis will differ depending on whether s is an integer or not, we will consider these two cases separately. Then we 
will allow the acceptable solutions to evolve up to the Cauchy horizon, and then on to the SFSH. 



1. Threshold, s . 



We can use the method of Frobenius to describe solutions near regular singular points, and we present these using 
the P-symbol notation ,23;]. For the first four columns in the P-symbol, the first row's entry denotes the location of 
the regular singular point, the second and third rows' denote the leading order exponents of two infinite power series 
solutions at that point. If these exponents do not differ by an integer, then these series are two linearly independent 
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solutions for D near that point; if the exponents do differ by an integer, then a logarithmic term may be introduced. 
Thus, for A < 1/8, 



2 (' 







V1-8A 







V1-8A 



(5.8) 



As this is a second order equation coming from a fourth order system, there must be two solutions with D — 0. We 
find these by setting D = in the system, which simplifies the equations greatly. Thus we can find the set of exact 
solutions corresponding to D = 0, which we call solutions /// and IV and are valid everywhere and irrespective of 
whether s is an integer or not: 





III 


A = 




D = 





K = 





G = 





IV 

is - l)kox'-^ - (£+^kox-'-i 


kox'-'^ 
0, 



(5.9) 



where fco,5o are arbitrary constants. The full set of solutions are found using system methods, as this approach has 
the benefit of solving for all four solutions at once, with the required accuracy found by simply taking further terms 
in the series. These system methods are described in Appendix C, and we briefly outline their application here: 
The system of perturbation equations can be written in the form 



Y' = ^(j+f:ArnxA 
\ m=l / 



(5.10) 



where Y = {A, D, K,G)^ , and J 0. Thus x = is an irregular singular point. We find that J has eigenvalue 
multiplicity four. J cannot be diagonalised and we use Theorem IC . 41 given in Appendix C to remove off-diagonal terms, 
which effectively reduces the order of the singularity to a regular singular point. Now the leading order coefficient 
matrix has eigenvalues 0,s,s,s — 2, and so we apply Theorem IC.2I twice to reduce the eigenvalues to and s — 2 
multiplicity three. Finally we apply Theorem IC . II to obtain the following. 

To leading order, we find a full set of linearly independent solutions with asymptotic behaviour (as x — + 0) 





/ 


II 


/// 


IV 


A = 


0(1) 


0{x'-^) 


Oix") 


0(x^-2 


D = 


0(1) 


0(x''+2) 








K = 


0(1) 


0{x'+^) 





0(x"-i 


G = 


0(1) 


0(x^+2) 


0(2;") 






Solutions / and // correspond with the Minkowski solutions (matching across a; = is dealt with in the next 
subsection), and III and IV are the D — Q solutions given in (|5.9(l . 

Taking A and D to be a linear combination of these solutions, we calculate the leading order of the scalars near 
x = as 

5*0^O(1) + O(:e^+2), 5-^i^0{l) + 0{x'-^)+0{x'), (5.12) 

(for brevity we have left out the constants of combination). Thus for our master function 5P-i to be finite on the 
threshold due to the Vaidya solutions, we find we must maintain the same constraint as that coming from the flat 
space solutions, that is 

Re{s) > 2, s ^ Z. (5.13) 



2. Threshold, s G Z. 

When s is an integer the system methods break down for the following reason: the eigenvalues of the leading order 
coefficient matrix of the regular singular point at x = are and s — 2 multiplicity three. These differ by an integer, 
and thus they must be repeatedly reduced until they are equal. However, each time we reduce an eigenvalue we 
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must diagonalise the leading order coefficient matrix, which prevents us from simply reducing the eigenvalue (the 
unspecified number) s — 2 times. 

Instead, we use the ordinary differential equation for D to write the solutions near a; = as, using the fact that 
s + 2 > (since the flat space solutions constrained s > 2), 

oo ( oc oo ^ 

D = d^Y^ + 4 i yfclna; ^ A„a;™+'*+2 + ^ B^x"' \ , (5.14) 

m— V m— m— ) 

with o?5, constants and k possibly zero. 

We use this solution as an inhomogeneous term to solve for the other variables by integration, and we find 

^ ^ ra + 2 \ ^ m- s ^ m + 2 \ m + 2/ |' ^ ^ 

V , ^ (-2(m + . + 2) + x(l-.-2A)) +dJ Vi?>n:r'" ^ ^ 

- E (^rli - (- - ;;r^) - - ;^)) } • 

m— ^ ^ ^ ^ ) 

{s + P + l-l)x-^D + 2\G+{s^ -l)x-^K - sK'. (5.17) 

Since both 13 and K have 0(1) terms, we see there is an A solution which diverges on the threshold like x~^. This 
divergent term cannot be switched off, for the following reason: 

On the axis there were two solutions for D, which we denoted 

^D,^x'-''\ ^D2^x'+K (5.18) 

The scalars 5^o,4 due to the first solution went like r'~^, whereas the second solution gave (55'o,4 ^ r~'~^. Thus we 
needed to switch off the divergent term in the general solution, -D|2;=-oo = j^Di + d2 A-D2, by setting c?2 = 0- 
Now when this solution was allowed evolve to the threshold, D\x=q — c?3 tDi + c?4 T-D2, the constants d^, c?4 7^ were 
fixed [see (|4.13|) ]. To match across x = 0, we must have (since we are using a global coordinate system) 

D^'\,=o = D^l^^o, (5.19) 

where M denotes solutions coming from Minkowski space-time and V denotes solutions coming from Vaidya space- 
time. Thus we require 

lim(d3 jD^' + d4 T-Df ) = lim(d5 ^.D^ + 4 t^*^)- (5.20) 

From H4.12|l . (|4.17|l . we see the solutions from Minkowski space-time are 0(1), 0(a;*+^). When s e Z > 2, we see from 
(|5.14|) the solutions from Vaidya space-time are also 0(1), 0(a;*+^), and thus to match continuously across a; = we 
cannot switch off the 0(1) Z3-solution. 

Thus when we calculate A, and hence S'i'4, there will be divergence as a; | due to this solution. This does not 
happen when s ^ Z, since there is no divergent ^-solution when Re{s) > 2. 

There were four constraints for initial data: (i) finite fiux on the axis, (ii) finite fiux on the threshold when approached 
from flat space, (iii) finite flux on the threshold when approached from Vaidya space-time, and (iv) continuous matching 
across x = 0. The most general class of perturbations which satisfies all of these conditions are those with 

Re{s) > 2, s^Z. (5.21) 



3. Cauchy Horizon. 



When A < 1/8, the Cauchy horizon is a regular singular point of the system given in Appendix B. Its leading order 
coefficient matrix has eigenvalues multiplicity three and 
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When a is not an integer, we can use the system methods outhned in Appendix C. Applying Theorem lC.il we find 
solutions with asymptotic behaviour 



/ 
0(1) 
D= 0(1) 
K= 0{w) 
G= 0{w) 

as w ^ where w ^ x — Xc (for consistency see H5.8II 
Now we make an important observation: since 



II 


/// 


IV 


0{w'') 


0(1) 


0(1) 


Olw") 













0(1) 


0(w'"+i) 


0(1) 






(5.23) 



< VT^8A < 1 



for < A < 1/8, therefore 



VI -8A 

and thus Re{a) > for Re{s) > 2. Alternatively, we can say 

cr = s - 2 + 0(A) 



> 2(2^-4), 



(5.24) 
(5.25) 
(5.26) 



where each coefficient of A" is positive, and thus again Re{a) > for Re{s) > 2. Thus each solution for A and D as 
given in (|5.23|l is at most 0(1) near x = Xc', all the solutions for A and D which are series beginning with w, w'^ or 
■u;<^+i will decrease to zero as we approach the Cauchy horizon. 
Since 



(5*0 = 



Ax - 1 ' 



r''-^{2A+{l - \x)D) 



(5.27) 



and A and D near the Cauchy horizon are a linear combination of 0(1) solutions, the scalars 54'o,4 representing the 
flux of the perturbation, and hence the scalar SP-i, will be finite on the Cauchy horizon x = Xc- Thus when cr ^ Z, 
the Cauchy horizon is stable under metric and matter perturbations. 



However, for each value of the parameter A < 1/8, there will be a mode number s such that cr G Z, and thus we 
must also consider this case. From (|5.8|l . we see a general solution for D near w = x — Xc = can be written as 



00 

E 



D = y ^ ArrrW 



d4 < /clnu; ^ ^mw"'^ 



m— 



rn— 



(5.28) 



where d^jd^ are constants and k can be zero. Since we are considering A > (A = being vacuum space-time) and 
Re{s) > 2, we have cr > 1 if cr G Z. Now we use this solution for D as an inhomogeneous term to integrate the 
perturbation equations. Near w — 0, we find a 4-parameter set of solutions: 



K 
G 
A 



9q 



1 — XXc 
1 

XciXXc - 1) 
1 



wD'dw — x^ ^{xc + 3 — n) 
Ddw, 



Ddw, 



D + 2XG + ^x-^[Xxc - Xc{l'^ + l-2)+6]K~ \{\xl + A)K' , 



where fco,(?o are constants, a prime denotes diff'erentiation w.r.t. w, and 



(5.29) 



Ddw = rfa E 

r 



°° AmW"'+''+^ I, 1 



m— 

00 



wD'dw ~ da ■ 

m— 



TO + cr + 1 

^™(TO + cr)u;"+'^+l 
TO + cr + 1 



m=0 



TO + cr + 1 



Inw — 



TO + cr + 1 



^ BmW 



m+1 



?TI = 



TO 



+ d4 < fc X! 



AmW 



,m+cr+l 



TO + cr + 1 



1 + (to + cr) In w — 



TO + cr 
TO + cr + 1 



C30 

E 

rn—0 



Bm mw 



rn+l 



1 



Since a > 1 and limu,^o w°' \nw ~ 0, we see all of these variables A, D, K, G, and thus the scalars (5*o,4 and SP^i, 
are again finite in the limit w 0. 

Thus the full set of perturbations which are finite on the axis and on the threshold Af will evolve up to the Cauchy 
horizon and beyond without their flux diverging. Therefore in the case of self-similar null dust there is a naked 
singularity whose Cauchy horizon is stable under metric and matter perturbations. 
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4- Second Future Similarity Horizon. 

Now something interesting happens when we aUow the solution to evolve past the Cauchy horizon and on to the 
next singular surface, the SFSH given hy ^ + ^/\ — 8A). The first scalar depends only on D, 

'5*0 = ^^^D, (5.30) 
Ax — 1 

and the solutions for D near x = Xe can be found directly fi'om H5.8|l as 

D^d^Y. ^™(^ - ^-)" + X] ^™(^ - ^eT^'^ (5-31) 

m— m— 

where 

^ = ^ f s - 4 - V (5.32) 



Since < \/l — 8A < 1, we see <j will always be negative for Re{s) > 2. Thus there is a class of solutions which are 
finite on the axis, finite on the threshold Af, finite on the Cauchy horizon, and then finally diverge on the SFSH. We 
emphasize that this instability is due to a; = Xg being a similarity horizon of the space-time, and not an event horizon. 



5.B. 1 = 1 mode, Vaidya. 



In this sector we can only define partially gauge invariant objects. As in Minkowski space-time, Section III.B, the 
metric perturbation objects are gauge sensitive to £,^dx^ = ^AYdx^ + £^Y,a dx"^ as 

kAB ^ [r'(eA'),A]|s+ [r2(CA'),B]|^ (5.33) 

k fc + 2e/r2 + 2wV(e/r2),^. (5.34) 

We examine the case where the perturbed space-time is that of null dust, and thus the bare matter perturbations are 
as given in (|5.2() . 

^^AB = L^J,-' ^..'^'-A.a.?^ , At^ = _8,-r , Ati = At' = 0. (5.35) 



Symm x Sp- y V ^ , 

Since pA = hA in this sector, we can transform into the equivalent of the Regge- Wheeler gauge by choosing 

U^hA^r'iar%A, (5.36) 

the benefit being that in this gauge pA = and therefore Tab — AtAB etc. Thus we can express the right hand side 
of the perturbation equations given in Appendix A in terms of 6p, T. 

Further transformations maintain this condition provided ^a = ^^^(^/''^)|A- Importantly, this fixes £^a while keeping 
^ completely free. 

Tab and Ta are not gauge invariant, they are sensitive to gauge transformations as 

Tab - Tab - iABic^'ieA')"^ + tcB(r2(eA')l^)|A + icA{rH^ynB, 

Ta-Ta = -?ABr2(e/r2)l^. (5.37) 

The vanishing of and is gauge invariant. 

Now we look at the perturbation equations in {x,r) coordinates. As in the I > 2 sector we perform a Mellin 
transformation of the equations, which is equivalent to parameterizing the perturbation components as in H5.3|l . 

kAB = ( ""^l^f (if r'-?c[t) ) ' ^ ^ ^'^'Ksixh r = r^Gsix), dp = r^-Hl,{x). (5.38) 



Again, dropping the subscript s, we define the new variable D = B — xA 

„2^ ,A 



We exploit the gauge freedom and transform into a gauge in which k^"^ = 0, by choosing ^ such that 



[rH^/r%Ar = -kA''. (5.39) 
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Then we are allowed make further transformations which preserve RW gauge and fc^"* = provided 

VHi/r^A]^^ = 0, 
U + r\i/r^)\A - 0. (5.40) 

Thus we recover the perturbation equation (|A.ldp . which was not valid in the I = 1 sector, as a gauge choice. As 
before the set of perturbation equations in this gauge reduces to one constraint equation defining 5p, a first order 
system in A, D, K and G, and we can decouple a second order ordinary differential equation for D. 



There is some gauge freedom in ^ left, and we will use this to gauge away k. Let us formalise this: let ^ satisfy 
gauge conditions Li^ = 0; let k satisfy perturbation field equation L2A; — 0; and let k transform as A; — > fc + La^. 
Since ^2X3^ = subject to Li^ = 0, there is a gauge in which k = Q. 

In the coordinate system (x,r), where we perform a Mellin transform over r such that ^ = r''~^^s, the above equations 
are 

LiC = -(-1 + s) (-.tA + s (-1 + xA)) is(.x) + {-2 + 2x-^x^ \ + 2s {\- x + x^ X)) ^'^{x) 
+x {-2 + x^x^\) £,'^{x), 

L2(. = C'(^)+Pi(^)C(^)+9i(^)eU^)+^i(^)6(^), 

= {s + xX~ sxX) £,s{x) + {l-x + x^ \) £,'^{x), (5.41) 

where 



qi{x) = 



6 + 5 (-2 + 3;) x-Ax (-3 + a; + 3a;^) A + a;^ (24 + 72:) + (-l + a;A) (4 + 3a; (-1 + 2:A)) 
X (l + n-x + x (3 - n + x) A) (2 + a; (-1 + a; A)) 
n (2 + z (1 - 16A + Z (-3 + A (le + T (6 - (17 + 3a:) A))))) 
a:: (1 + n - a; + a:: (3 - n + a::) A) (2 + X (-1 + a:: A)) 
(2 - n) (-2n(l +n) + (1 + n) (-1 +3n)a;+ (1 - 3n) a:^) - 2a; (3 + n (-1 + (-3 + n) n) - a; - (-3 + n) n (-4 + 3n) a;) 



+ 



+ 



r\{x) = 



a;2 (1 + n - a; + a: (3 - n + a;) A) (2 + x (-1 + a; A)) 
(-2a; (5 + 3 (-3 + n) n) a;^) A+(-l + n) (-30-8a; + n (19-3n + 3a:)) A^ 
a;2 (1 + n - a; + a; (3 - n + a;) A) (2 + a; (-1 + a; A)) ' 
(-1 + ri) ((2 - n) (1 + n) (n - a;) + 2 (3 + a; + a; (5 + a;) + 2 n (l + a;^) - (1 + x (4 + x))) A) 
x2 (1 + n - x + X (3 - n + x) A) (2 + x (-1 + x A)) 
(-1 + n) ((-3 + n) n (-4 + n - x) x^ A^) 



x2 (1 + n - X + X (3 - n + x) A) (2 + x (-1 + x A)) ' 



A direct consequence of fc = is that D = 0. Thus we are in a gauge in which pA — kj^^ — k — D — 0, and to 
remain in this gauge we are allowed further transformations provided 



2^/r' + 2i''r'{^/r'),A = 0. (5.42) 

The only ^ which satisfies both these constraints is ^ = 0. Therefore there is no remaining gauge freedom and so the 
remaining perturbation variables are gauge invariant. 
Thus we have found a one parameter family of solutions, 



, _ { 2goXr^+'x^ 2g,Xr^x^+^ \ ^ _( x^s+^j^^^ 
\2gaXr'x'+^ 2gQXr'-^x'+'^ ) ' \ \ao{s+x)T'--'x' \go(s+x)r' 



3^s + l 



47r 



^^^j, Ti=r2 = /c = 0. (5.43) 

To match continuously with the empty Minkowski region preceding x = we require Re{s) > 1. These solutions will 
evolve without divergence through the rest of the space-time. 

Note the perturbations will vanish if go = 0, that is to say if we had considered only metric perturbations, and no 
matter perturbations, we would have returned an empty sector. Note also the sector is empty in the Minkowski limit 
0. 
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5.C. I — mode, Vaidya. 

When / = 0, then Y,a — and our perturbations are 

h f hAB \ ( ^tAB \ 

We will use the coordinates {v,r,9,(j)), where v is the null coordinate of the background. Thus the metric of the 
background is g^^dx^dx'^ = —(1 — ^)dv'^ + 2dv dr + r'^jabdx'^dx^. 

Since the matter tensor has the form H5.1() and has no angular dependence, we find At^ = 0. We can describe 
the ingoing radial null geodesic of the perturbed space-time as + S£f_i — — V^F where V = v + T{v,r) is the null 
coordinate of the perturbed space-time. Our unknowns therefore are Hab, K, Sp and T. 

We cannot construct gauge invariants in the / = sector and thus we exploit the remaining gauge freedom to set 
some variables equal to zero. The perturbation variables are gauge dependent as 

hAB hAB - {S,A\B + ^B\a), 

K ^ K-2i^U. 

^tAB ^tAB — iAB\C^'^ ^ icBi.'^\A ^ ^CaC^iB' (5.45) 

As before, we transform into a gauge in which K = /i^^ = 0. To do this we choose such that (in the (t;,r) 
coordinate system where 9 and 9' denote differentiation w.r.t. v and r respectively) 

e. + (l-4^)er = \rK, 



ir = hh^A- (5.46) 



2 
1 
2 

Then we are free to make further gauge transformations which preserve this condition provided 

+ (1 - ^)ir = 0, 

= 0. (5.47) 

Now we look at the field equations in this gauge (we will let h^ — A, hyr — hrv — B and h^r = C). The first is 

(7 = 0, (5.48) 

and thus C = C{r). When we perform a gauge transformation on this quantity subject to H5.45|l and (|5.47|l . wc find 
C — > C — 2^^, but since is an arbitrary function of r, this means we can choose a gauge in which C — Q (and thus 
B = Q since h^j^ = 0). Thus we have transformed to a gauge in which K = /i'^^ = B = C ~ 0, and to remain in 
this gauge we are allowed further gauge transformations of the form £^a = co(— (1 — ^)(5^ + (5J^), with cq an arbitrary 
constant. 

The remaining perturbation equations in this gauge are 

rA" + 2A' = 0, 
rA' +A + XT' = 0, 
^{1- ^){rA' + A) + lA~^t = 87T6p. (5.49) 

That + S£f^ must be null and geodesic gives F' = 0, and hence 

r = a{v), A=^, 5p^ -^{l3-2Xa). (5.50) 



Further gauge transformations give 



A A~2^co, 

At AS AtAB, (5.51) 



and thus these remaining perturbation quantities cannot be gauged away. 

What we have shown here is essentially a uniqueness result: all the above perturbations can be generated by a 
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perturbation in the mass function and the null vector. The metric and matter tensors for spherically symmetric null 
dust are given by 

gf.^dxf^dx" = - (^1 - _^ 2dvdr + r'^dn^, 

t^.^dx'^dx'' = (5.52) 

where before perturbation m{v) = Xv and = —d^v, and after a perturbation miv) = Aw + I3{v) — 2Aco and 
= —dfi{v + a{v)). The terms a,/?, cq are arbitrary and therefore there is no reason to suspect divergence on the 
Cauchy horizon or elsewhere. 



6. RESUMMATION 

In this section, we review some properties of the Mellin transform (see e.g. |2^, [S^) and present a plausibility 
argument that finitcncss of the modes (as seen above at the Cauchy horizon) implies finiteness of the full perturbation 
found by resumming the modes. This argument is based on a study of the Klein-Gordon equation (or wave equation) 
in Vaidya space-time, and its correspondence to the wave equation in Minkowski space-time. 

We recall that the line element of self-similar Vaidya space-time may be written as 

ds^ = r2(-l + Xx)dx'^ + 2r(l - x + \x'^)dx dr + x{2 ~x + \x'^)dr'^ + r'^dVL^, (6.1) 

with the coordinate ranges r e [Q,oo),x e (— oo,cx)). We note that Minkowski space-time corresponds to taking 
A = 0. Then the wave equation (or more accurately, the PDE satisfied by the Z**^ multipole moment (f) = ipi oi the 
Klein-Gordon field) is 

- x{2 - x + Xx^)(j)^xx + 2(1 - a; + Ax^)r0,„ + (1 - Xx)r^(p_rr - Xx^cf)^^ + (2 - Xx)r(t)^r - l{l + l)4> = 0. (6.2) 

The Mellin transform of this equation yields a parametrised ODE. The first step is to define the Mellin transform of 
the field (note that the arrow indicates the variable with respect to which the Mellin transform is taken) : 

t-i 



P{x;t) = M[^{x,r)]{r ~^t) -.^ (j){x,ry-' dr, (6.3) 

"'0 

where t £ C and satisfies ri < Re{t) < T2, the constants ti^2 being defined by the condition that the integral converges 
for these values of t. The field (l){x, r) is recovered via the inverse Mellin transform: 

r+zoo 

{x;t)r-*- dt, (6.4) 



{x, r) = M-^[P{x; t)]{t ^ r) = — P 



where the inversion contour is vertical and lies in the strip ri < Re(t) < T2. We should note here that as the 
resummation is over terms of the form P(x;t)r~^ , the role played by t in the present section corresponds to the role 
played by — s in Sections 2-5. See in particular the last paragraph of Section 2. 
Integrating by parts immediately yields the following results: 



Lemma 6.1 



provided 



Lemma 6.2 



M[r(l)^rix,r)]{r t) = -tAl[0(x,r)](r t) = -tP{x;t), (6.5) 



lim r*(/)(x, r) = lim r*(j){x,r) = 0. (6-6) 

r — *0 r — ^CJO 



M[r'^ (I), rr{x, r)](r ^ t) = t{t + l)M[(j){x, r)]{r t) = t{t + l)P{x; t), (6.7) 
provided i)6'.6|) holds and 

lim r'+V.r(a;, r) = lim r*+ V r(a;, r) = 0. (6.8) 

r—*0 ' T — >oo ' 
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Assuming that these conditions hold for values of t on a strip of the complex plane, we can take the Mellin transform 
of iniSl to obtain the ODE 

-x{2-x + Xx^)P" - (2i(l - x) + (1 + 2t)\x^)P' + {t^ -t- l{l + 1) - f\x)P = 0, (6.9) 

where the prime here and throughout refers to differentiation with respect to argument. This equation has analytic 
coefficients (and hence analytic solutions) everywhere except at infinity and at the singular points x = 0, XcXg, the 
roots (in increasing order) of x{2 — x + \x^) — 0. Note that in Minkowski space-time, there is no third singular point 
Xe- These are all regular singular points of the equation, and thus the standard Frobenius theory can be used to study 
the global behaviour of solutions of H6.9|l (see for example (2^], or any textbook on linear differential equations). We 
recall that x — xq ^constant is a space-like hypersurface for xq £ (0, Xc) and that x — 0,Xc are null hypersurfaces. 

For Minkowski space-time, A = and H6.9|l is a hypergeometric differential equation, and we can give an essentially 
complete account of the problem at hand. We proceed to do so in order to clarify the nature of this problem and our 
putative solution. So let A = in (|6.9(l and for convenience of comparison with the standard text by Bateman [2^ . 
let X — 2z (all results quoted below are taken from this reference). Then (|6.9|) is the hypergeometric equation 

z(l - z)u" + (c-{a + b+ l)z)u' - abu ^ 0, (6.10) 

where u{z) = P{x), a = t + l,b = t — I — 1 and c = t. We will assume that i ^ Z. This can be assumed without loss 
of generality by a deformation of the inversion contour. 

The past and future null cones of the origin then correspond to the singular points z = 0, 1 respectively of this 
equation. We encounter here a slight difficulty. Being a null hypersurface, x — z — cannot be an initial data surface 
for the equation (|6.2fl . We expect that this will translate into z — failing to be a 'good' initial point for the ODE 
()6.1Q(I . However, our overall aim is to argue that finiteness of the field on x = along with finiteness of the modes 
at the future null cone (Cauchy horizon) is sufficient to imply finiteness of the field at the future null cone. We 
can connect these two by determining their respective connections to Cauchy data on a space-like hypersurface, for 
example, x = l(z = 1/2). So consider the Cauchy data 

We assume that a, (3 satisfy unspecified differentiability and integrability conditions that, in particular, allow us to 
calculate the Mellin transforms 

a{t) = M[a{r)]{r t), b{t) = M[P{r)][r ^ t) 

on some strip ti < Re(i) < T2 of the complex plane (the factor 1/2 is given for later convenience). These then yield 
the initial data for the Mellin transform P{x\ t) oi 4> &% x = 1, i.e. for u{z) at z = 1/2: 

n{l) = a{t), u'i^)=bit). 

We emphasise that our assumptions on a, /? imply the existence of the inverse Mellin transform of a, b taken over a 
vertical contour in ri < Re(i) < T2. 

The next step is to determine the solution for m at 2 = and at z = 1 in terms of a{t) and b{t). We use the following 
pairs of linearly independent solutions at these two points (again following the notation of Bateman ^^). At z = 
we use 

ui{z;t) = (1- z)^~'^F {-1,1 + 1; t;z), U5{z;t) = z^'*' F{1 + 1,-1,2 - t; z), 

and at z = 1 we use 

U2{z;t) = zJ-^FH + 1, -l;t; 1 - z), UQ{z;t) = (1 - z)'^-'^F{-l,l + 1;2 - t; 1 - z), 

where 

F{a, b; c; z) = ,F,{a, b; c; z) = M^z", 

with 

(a)o = 1, (a)n = a{a + 1) • • • (a + ?i - 1), n = 1, 2, 3, . . . 



20 



is the standard hypergeometric function. Note then that all the hypergeometric functions present are in fact polyno- 
mials, and so M5 and ug are finite sums of powers of z and 1 — z respectively. U\ and 1*2 are infinite series in z and 
\ — z respectively, both with radius of convergence at least 1. We will concentrate on the solution at z = 0. The 
general solution of H().10|l at z = is 

u(z;i) = ci(t)Mi(z;t) + c^(t)u^{z\t), (6-11) 

giving 

u'(z;i) =Ci(tK(z;t)+C5(t)u'5(z;i). (6.12) 

Note that existence of u^\z=ii and U5|z=o requires Re(<) < 0: ruhng out the integer case, we take Re(t) < 0. As 
z = 1/2 is within the radius of convergence of the series solutions, we have 

ci(i)wi(^;0 +C5(t)u2(i;t) = a(t), 

ci(iK(^;0+C5(t)4(i;0 = hify. 

Solving for ci(i),C2(t) gives 

CiW = ^(a(t)u:,(i;i)-6(iK(^;i)), 

C5(i) = ^(-a(tK(^;i)+Ki)«i(^;i)), 

where we have used Abel's formula for the Wronskian which arises as a determinant in solving for ci,C5. We can 
obtain similar expressions for C2(t\c^{i)i the coefficients of M2(z; t), U6(z; i) in the general solution for u at z = 1. 

Our next step is to address the following question. Do the conditions on a{t),b{t) (i.e. existence of the inverse 
Mellin transform in the strip ri < Re(t) < T2) imply the existence of the inverse Mellin transforms of 

P{0;t) = c,{t)ui{Q;t)+C5{t)u5{0,t), 
P'iO;t) = i(ci(tK(0;t)+C5(t)<(0,i)), 

and of the corresponding expressions at a; = 1? By linearity, this is true if and only if it is separately true for a and 
for b. We consider only the case a ^ 0, b — Q: the reverse case is similar. Then noting that ui{0;t) = 1 and that 
1/5(0; i) = u'^{0;t) = (since Re{t) < 0), we have 

to{t) := P{0;t)^^u',{^;t)a{t) 

= £^^^'(t-l)(t-2)...(^-fc)«(^) 



where 



CM = (-l)''|j(' + l)fe(-Ofe, A: = 0,1, 



The rational function of t can be written as the sum of inverse linear terms using 

1 (-l)'^-J' 1 



(t-l)(t-2)...(i-fc) ^^(fc-j)!(j-l)!i-j- 

We now point out two properties of the Mellin transform that will allow us to perform the required inversion. 
Lemma 6.3 If 

Git)^MMr)]ir^t), 

then 

for < fc e M, provided all relevant integrals converge. 
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Proof: This is immediate from the formula (|6.4() for the inverse Melhn transform. 
Lemma 6.4 If 

then 

for fc S M, provided all relevant integrals converge. 

Proof: This follows from Lemma 6.1. 

Using these results, we could give a closed form expression for the inverse Mellin transform of to{t) in terms of 
integrals of a(r), the inverse Mellin transform of a{t). These will have the structure of a (terminating) series, whose 
A:*'* term consists of a sum of k terms, the j*'* of which has the form 

^^^1^ I ' ^{y)y'^^dy (6.13) 

for some constants dj^k- Note that the r/4 arises from the term 2~'^* in to(^) sjid by using Lemma 6.3. Since there is 
only a finite number of terms present, the only questions regarding convergence relate to the existence of the integrals 
(|6.13|) . This existence follows from the conditions laid down on the initial data function a(r). In exactly the same 
way, we can deduce the existence of the inverse Mellin transform of 

ti{t) = P{2-t) = u{l-t) = '^——u'^{l-t)a{t), 

and of the terms corresponding to io(i) and ti(t) that arise by taking a{t) = 0, hit) ^ 0. 

To conclude the discussion for Minkowski space-time, we take a slightly different point of view, where we treat 
u{Q;t) — ci(t)ui(0; t) + c^{t)u^{Q] t) and u'(0; t) — ci{t)u'i{Q] t) + C5(i)u5(0; t) as the fundamental data for the problem. 
Choosing ci{t), c^it) so that the inverse Mellin transforms of these terms exist on a strip of the complex plane implies 
(by the work done above) that the inverses of u and u' also exist on the space-like hypersurface z = 1/2. Then 
verifying that the restriction on t implies finiteness of the basis solutions U2, ue at z = 1, we can finally conclude that 
the inverse Mellin transforms of P and P' exist at the future null cone z — 1. 

The thrust of this argument is that the 'scattering coefficients' Ci{t) do not affect resummability of the solution, 
and so that finiteness of the basis solutions at z = 1 is a necessary and also sufficient condition for finiteness of 

the field (f> at the future null cone. 

Apart from actually verifying these assertions rigorously, the last step in our argument involves showing that the 
scenario summarized in the two preceding paragraphs carries over to Vaidya space-time. The crucial step is to consider 
the dependence of pairs of fundamental solutions of the ODE at the two crucial singular points: {Pi tP^ ) at the 
past null cone and (pf ,^2') at the Cauchy horizon (future null cone). We consider first solutions at x = 0. 

The standard form of H6.9|l at a:: = is 

x^P" + xq{x; t)P' + r{x\ t)P = 0, 

where 

2t - 2tx + + 2t)Xx^ 
2-X + A.t2 ' 
{t^-t- l{l + 1) - t^\x)x 
2- x + Xx"^ 

Note that these are analytic at a; = and can be written in the form 

00 oc 

q{x;t) = ^gfe(^)a;^ r{x;t) = ^rfe(^)a;^ 

fe=0 k=Q 

where convergence of the series is guaranteed for x G [OjXc). The indicial equation is 



q{x;t) = 
r{x;t) = 



I{v) := v{v - 1) + qav + r^ ^ v{v - I + t) ^ Q, 
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with solutions i/ = 0, 1 — f ^ Z. Thus we have the linearly independent solutions 

oo 

p^{x;t) = '^ak(,t)x'', ao = 1, 

oo 

p^{x;t) = \x\^-'Y,h{t)x'', bo = l. 

In each case, the radius of convergence of the series is at least Xc- Thus this representation of a fundamental set of 
solutions is valid for all x € [0,Xc)- The recurrence relations for the ak{t) are 

^ fe-1 

Mi) = -jm^Ulk-jit) + rk-j{t))aj{t), k>l, 

and the recurrence relations for the bk are 

^ fe-1 

Wc wish to determine the f— dependence of these coefficients. Noting that the coefficients qk{t) and rk{t) are respec- 
tively linear and quadratic in t and that I{v) is linear in it is straightforward to prove that 

where Vk{t) is used to represent an arbitrary polynomial of degree k which may be different in different formulae and 
within a single formula. Then division in the ring of polynomials yields 

with a similar result for 6^. From here it is clear that we can write more explicitly 



au = Vkit) + ^k-ijt) ^ y jr^ 



t{t+l)---{t + k-l) t^o^-^' 

= + (2 - m + 1 - ^) = ^'^'^ \§ Jrf^' 

where the a-j. and bj. are independent of t. 

Before proceeding, we note another property of the Mellin transform that will allow us to deal with the inversion 
of terms arising from the presence of the Pk{t) in the coefffcients ak{t). 

Lemma 6.5 // 

Git)^MMr)]ir^t), 

then, provided the relevant integrals converge, 

s'^Git) = M{D''-f{r)]{r ^ t) 

where the differential operator is defined by D = —r-^. 

Proof: This follows by induction and by using Lemma 6.1. 
There is an immediate corollary: 

Corollary 6.1 // 

G{t)=MMr)]{r^t), 
and p{t) — po + Pit + ■ ■ ■ +Pnt'^ is a polynomial in t, then 

p{t)G{t) = M[(po+PiD + • • • +p„i^")7(r)](r ^ t). 
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Note that the preceding analysis apphes equally well to the case A = 0, i.e. to the wave equation in Minkowski 
space-time. Thus the solutions Ui,u^ can equally well be written as combinations of We emphasize this as 

we now have solutions in Minkowski space-time that involve infinite series rather than polynomials. 

The solutions , i = 1,2 'scatter' to a naturally arising fundamental set of linearly independent solutions Pi ,P2 
defined at a; = a;c (where = 2 in Minkowski space-time). Indeed by writing (|6.9I) in standard form aX x = Xc and 
determining the roots of the indicial equation thereat, we can write explicitly 



Pi 



\x-t) = (6.14) 



n=0 



p'^(,x-t) = \x~Xc\'''Y.^k{t){x-x,T, (6.15) 



n=0 



where 



and both series converge in some neighbourhood of a; = Xc- We note that Re(;/2) > Re(l — i), and so the earlier 
restriction Re(i) < implies that Re(i/2) > 1, and so both solutions = 1,2 are finite at the Cauchy horizon. 
(This applies quite generally in self-similar collapse to a naked singularity when the dominant energy condition holds: 
see 0.) As these solutions can be written as different C— linear combinations of p^ ,1 = 1,2, this implies that the 
series representations of pf' ,i = 1, 2 must both converge at x — Xc- This is of crucial importance to our argument, as 
we can now determine the field (pix, r) at the future null cone x = Xc hy carrying out (i.e. checking convergence of) 
the inverse Mellin transform of the sum of the convergent series 

Pixc] t) = a{t)p^{xc; t) + p^[x,- 1). (6.16) 

Here a(i), 6(t) are functions whose inverse Mellin transforms a(r), (3{r) exist for inversion contours lying in some strip 
of the complex plane, a, (3 constitute Cauchy data for 0, and satisfy unspecified differentiability and integrability 
conditions. 

As we have seen, this inverse Mellin transform will involve an infinite series of terms of the form 



Ctk 



(r) = Vk{D)a{r) + ^ fe,r=- f ^ dy, 

^ h y ' 



where D is the differential operator introduced above. There will be similar terms arising from the contribution 
by P{r), the inverse Mellin transform of b{t). Existence of these individual terms can be guaranteed by imposing 
differentiability and integrability conditions on the initial data functions a(r),/?(r). (Assuming differentiability of 
arbitrarily high order should not be a significant constraint here as linearity would allow us to work in distributions or 
to use a density argument to generalise from analytic functions to more interesting spaces. We also note that the sign 
of the Ci in the above will not be of particular relevance, as for either case of this sign, either the multiplicative 
pre-factor or the divisor in the integrand will have a mollifying effect.) Hence our problem boils down to deducing 
convergence of series of the form X^fcLo '^ki'f)- Here is another crucial point: our analysis of the problem in Minkowski 
space-time from a slightly different point of view guarantees that such series must converge in the case A = 0. 

Now the coefficients of these series are generated by the coefficients of the series in ()6.16(l . The only difference 
between Minkowski space-time and Vaidya space-time is the value of A (A = and A £ (0, 1/8) respectively). For 
A = 0, the coefficients of the series of complex numbers (|6.16f) guarantees convergence in C: these coefficients generate 
coefficients of an infinite series + dkPkir)) in a certain unspecified function space - call it JF - that, as we 

have argued, guarantee convergence in this space. Convergence of this series depends only on the coefficients inherited 
from 1^6.16]) . To conclude our argument, we maintain that this pattern is repeated when A > 0: the coefficients of 
the C— convergent series (|6.16|l generate coefficients of series of functions in a certain unspecified function space that 
must then also converge in this function space. 

We conclude this section by summarising the argument. Restricting the values of t to allow only a finite fiux at the 
regular axis and at the past null cone is sufficient in the present case to guarantee finiteness of the Mellin transform 
P{x; t) of the field (j) at the future null cone. To determine if (j) itself is finite at the future null cone, we must 
calculate the inverse Mellin transform of P. The analytic form of the line element of Vaidya space-time allows us to 
determine the general form of this inverse: it involves an infinite series of finite sums of derivatives and integrals of 
the initial data for </>. As these finite sums converge, convergence of the full series of functions in depends only on 
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the constant coefficients of the series: these coefficients are generated by the coefficients in the convergent C— series 
representation of P{x;t) at a; = Xc- In Minkwoski space-time, the coefficients {Pk{0)}'kLo of the series P{xc\t) of 
complex numbers produces a convergent series representation for 4'\x=xc i^i ^- We claim that in Vaidya space-time, 
the same will happen: the coefficients {pfc(A)}^Q of the series P{xc',t) of complex numbers produces a convergent 
series representation for (j}\x=xc in ^■ 

This last paragraph constitutes our argument that in order to determine finiteness of the field at the Cauchy horizon, 
it is sufficient to check finiteness of the modes thereat. There clearly remains a good deal to prove, but the analysis 
above suggests a way of doing this: the main thing that requires checking is the convergence of the finite sums ak{r) 
and of the overall series J^'^k- Finally we note that although this argument has been presented for only the wave 
equation, it should generalise to any set of linear equations, as for example the perturbation equations considered in 
this paper. 



7. CONCLUSIONS. 



We have considered metric and matter perturbations of all angular modes falling on the Cauchy horizon formed 
by the naked singularity arising from the collapse of a self-similar null dust. There is no class of perturbation which 
satisfies the initial conditions and gives a divergent flux on the Cauchy horizon, thus there is no "blue sheet" instability 
as is seen in, for example, the Reissner-Nordstrom naked singularity. (We have shown rigorously that the modes are 
finite at the Cauchy horizon, and given a plausibility argument that the full perturbation itself, found by resumming 
over the modes, is finite.) Interestingly, the second future similarity horizon of the self-similar null dust space-time is 
unstable for perturbations with multipole modes I > 2. 

The question of uniqueness of the perturbation to the future of the Cauchy horizon has not been addressed, but this 
should not affect the divergence encountered at the SFSH. If we consider the wave equation as studied in Section 6 as 
a paradigm, we note that the solution (|6.15() and its derivative vanish at the Cauchy horizon for the allowed range of 
t (Re{t) < 0). Thus an arbitrary additional amount of P2 could be added to the solution for x > Xc while preserving 
continuity and differentiability. However in Minkowski space-time a selection procedure tells us the correct addition 
to make, based on reflection through the (still) regular axis. It is not clear that one can apply a similar argument in 
Vaidya space-time, where one encounters a singularity at r = to the future of the future null cone (Cauchy horizon): 
indeed this is the central problem caused by a naked singularity, and the very issue that the Cosmic Censor seeks 
to render irrelevant. Whether establishing uniqueness is possible or not, some amount of both independent solutions 
pf, i = 1,2 will persist in x > Xc, leading to divergence at x = Xe- 

It is worrying from the point of view of the Cosmic Censorship Hypothesis that this naked singularity persists. 
However, this worry is tempered by the fact that an instability is encountered at the SFSH of the space-time; the 
naked singularity only survives for a finite time. It is then of interest to consider the life-time of this naked singularity 
and what effects it might have on the space-time. It is also of interest to know how generic this short-lived stability 
of the naked singularity is. In the terminology of Carr and Gundlach |17| . this second future similarity horizon is a 
"splash" , whereas the other similarity horizons are "fans" . We speculate that finiteness of the flux of perturbations 
on a fan-type similarity horizon, and divergence on a splash-type similarity horizon, may be a general feature of 
self-similar spherically symmetric space-times; the issue is currently being studied. 

The results here also have some bearing on the issue of stability in critical collapse. In perfect fluid collapse, 
the critical space-times are continuously self-similar, i.e. admit a homothetic killing vector field (see e.g. |0|). By 
definition, the critical space-times posses a single (spherically symmetric) unstable mode for perturbations injected 
along the homothetic surfaces x — constant and followed up to the singularity. In our variables, this corresponds 
to a mode of the form r'^Q(x), with the limit r — > taken along x = constant, and with x < Xc- Then Re{s) < 
corresponds to instability. In the present paper, the limiting behaviour of the perturbation has been studied in the 
approach to the first and second future similarity horizons in Vaidya space-time. Carrying out a similar calculation for 
the critical fiuid space-times may reveal very different stability properties to that shown by critical collapse studies. 
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APPENDIX A: PERTURBATION EQUATIONS IN TERMS OF GAUGE INVARIANTS. 

We give here the full set of perturbation equations for the gauge invariants defined in ()3.14fl . Note not all equations 
apply for ^ = 0, 1 modes. 

2^^(fcAB|C - kcA\B - kcB\A + 2gABkcD^'') - (^-^^^^ + + Ga" + 2^) kAB 

f + \{Gc'' + G^) + n] k^^ 2gABv'^k^\, + g^B(2«^l^ + Av^v^ - G^^)kcD 



-~gAB ( 2k,c +6«^fc,c TT^fc) + 2(«aA:,b +iBk,A +k,A\B ) = -IGt^Tab, {I > 0) (A.la) 

-{kj^ +2i^k,c +G,-k) + (fcp^l^l^ + 25^fc^^l^ + 2(^^l^ + ^^5^)fccD) 

fc^^l^l^ + 7^fcc:^ - ^^fcc^) - -167rT\ (l > 0) (A.lb) 

k,A-kAc^^ + kc^^^-iAkc^ = -16TrTA, {I > 1) (A.lc) 

fc^^ = -^^WttT^. {I > 2) (A.ld) 



Here = and G^i, is the Einstein tensor of the background space-time. TZ is the Gaussian curvature of M^, 
the manifold spanned by the time and radial coordinates, and thus equals half the Ricci scalar of M^. 



APPENDIX B: FIRST ORDER SYSTEM OF PERTURBATION EQUATIONS. 

The perturbation equations give rise to the first order system Y' ~ M{x)Y where Y — (A, D, K, G)'^ and the 
coefficients of M are: 



Mil = 

Mi2 = 
Ml3 = 



Mi 



14 



4-4 x + 2 (l + X) -3 X x^ + \^ x^ + 2 {2-x + \ x^)+a [-4+2 x- \ x^ + \^ a;*) 
X (2s + Aa;2) {2-x + \x^) 

4 {-l+x)-2 X (-1+X + l (2+x (-1+1^ A)) (2+x (-1+x A) (2+x (-1+x A)))+i^ (2+x (-l+o: A)) (2+x (-1+x A) (2+x A))) 

x^ (-1+X A) (2 a+x^ a) (2+2! (-1+2^ A)) 
-(s (-1+X A) (4+x (-2+X (2+A (-2+3: (-5+3 x \))))))+x (8 \+x (-2+A (-12+3: (9+A (S+7 x (-2+x A)))))) 

x'^ (-1+X A) (2 s+x'^ a) (2+x (-1+X A)) 
(2+(-2+i + ;^) x) (2 + (-2+a:) x)-x^ (-6 + (9+2 I (1 + (-l + i:)-4a:) x) \ + x^ (3+(-2+i + i^) i:) A^ 

2:2 (2 s + a:2 a) (2+x (-l+o: A)) 
-2 (2+x (-1+X A)) + s a: (2-i: (2+A (4+3: (-5+3 x A)))) 

2:2 ^2 s + x^ a) (2 + i: (-1 + X A)) 
2 A (-4+2 s (2 + a: (-l + i: \))-x (-2+X A) (2 + a: (-l + o: A))) 
X (2 s+a:2 a) (2+x (-1 + x A)) 
_ 2j:(-1 + 3:A) 

(2s + i:2a) (2 + :r(-l + a:A)) (B.l) 
2 + 1 (l + l) (-2+a:)-((-l + s) (-6+4 s-x)+2 I (-l + 2:)+2 (-1 + x)) x A 
(-1+X A) (2 3+x'^ a) (2+x (-1+X A)) 
(6+2 s (-3+3-a;) + (3+!+!2) a;) \^ + (-2+s) x'^ \^+2 s (-2+3+2 A) 
^ (-1+X A) (2 s+x^ a) (2+x (-1+X A)) 

^ ^ (-2+2 3-(-2+i+i2)3;) (-1+3;A) ^ ^ _4^A(-l+a:A) 

(2s+a:2A) (2+a;(-l + a;A)) (2 s+a:^ a) (2 + a; (- 1 + a: A)) 

ji.f _ -2 » J- _ 2-2 s+l (l+l) (-2+x)-x (4-2 s+l (l+l) x) A 

^^31 - 27+^^ ^32 - .(-1+:.A) (2 s+x^ 

-2 + 2s2 + 3a:2A-a:(-2 + i + i2+a:A) 4 ;^ 

^33 = 277+Pn^ ^34 = 2s+a:^A 

M41 = M42 



M22 



M43 = M44 = J 
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APPENDIX C: METHODS FOR SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS WITH 

SINGULAR POINTS. 



For a first order system Y' = M{x)Y we define p as the least number such that the system can be written near 
a; = as 

and near x = oo as 

Y' = -xP-^ (j+Yl A„^x-'A Y. (C.2) 

V m=l J 

witli J 7^ and constant. 

We classify singular points and describe the solution as: 

p = 0, Regular point. The solutions of a system of differential equations are at least as regular as the coefficients 
of the system of differential equations. For the purposes of this paper, it is enough to know that solutions near 
a regular point are themselves regular, as none of the surfaces of interest are regular points. 

p = 1, Regular singular point. Also known as a simple singularity or a singularity of the first kind. Here we 
distinguish solutions depending on whether the eigenvalues of J given above differ by an integer or not. If they 
do not, we apply Theorem IC . II immediatelv. If they do, we reduce those eigenvalues individually until they are 
equal using Theorem IC. 21 and then apply Theorem lC.il (see ^28|). 

Theorem C.l In the system if6'.-?|) . if J has eigenvalues which do not differ by positive integers, then, in a disc 
around a; = not containing another singular point, fC'.ij) has a fundamental matrix $ of the form 

oo 

$ = Px\ where P{x) = ^ x™P™, Pq = I- (C.3) 

Theorem C.2 Let the distinct eigenvalues of J in \U.l\j he (disregarding multiplicity) /ii, . . ./i^, (k < N, where 
N is the order of the system). There is a matrix V{x) such that Y — VQ transforms f C. l\l into 

(C.4) 

where J has eigenvalues /ii — l,/i2, ■ • ■ iMfe- V is given by a diagonal matrix with entries 1, except where the 
eigenvalue(s) to be reduced occurs, where there is an x. 

p > 2, Irregular singular point. Also known as a non-simple singularity or a singularity of the second kind. Here 
we distinguish solutions depending on whether you can diagonalize J given in (|C.2|1 . If J has distinct eigenvalues, 
then J is diagonalizable and we apply Theorem lC.3l If J has multiple eigenvalues and J can only be reduced to 
Jordan normal form, then we apply Theorem lC.4l to remove off-diagonal terms (see "sS]). When the eigenvalues 
are repeated zeroes, this has the effect of reducing the order of the singularity, as happens in this paper at the 
threshold (see part V). 

There is a class of problems in between: sometimes a matrix has multiple eigenvalues and yet can still be 
diagonalized. In this case, there is a straightforward Theorem given by |30| if Ai = (as is sometimes the case 
when a high order equation is written as a first order system) . If not there is a very cumbersome solution given 
by 0. 

Theorem C.3 Let J have distinct eigenvalues //i, . . . , ^jln ■ Then !^C.1\) has a fundamental matrix 

oo 

•^ = Px^e", where P(a;) = ^ a;™P„, Po = /, (C.5) 

m— 

R is a diagonal matrix of complex constants, and H is a matrix polynomial (r — p ~ 2) 




H 



^^^_^Ha + yHi + ...+xH,, H, = dlag(^^x(^,...,^i^|^Y l^T = l^r (C-6) 
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Theorem C.4 For brevity's sake, we give the Theorem only for p = 2 as this is the case that arises in this 
paper. We transform J to it 's Jordan normal form J, and write the blocks of J as fil + pE, where E is the 
matrix with I's along it's super-diagonal and zeroes elsewhere. For each block of J, define the matrices 



D — diag{l, px, . . . , {px) ) 



B = 



I 1 1 1/2! ... 1/(A^- 1)! \ 
1 1 ... l/(A^-2)! 



V 



where N is the order of the system. Then the transformation Y = D ^BW gives the system 



(C.7) 



W 



pi + D'D-^ + B-^D I A„a;-™ ] D-^B 



W, 



(C.8) 



and the leading order coefficient matrix has had its off-diagonal terms removed. 
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